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We investigate how the anticipated neutral current rate from SNO will sharpen our 
understanding of the solar neutrino anomaly. Quantitative analyses are performed with 
representative values of this rate in the expected range of 0.8 — 1.2. This would provide 
a 5 — 10 cr signal for i/g transition into a state containing an active neutrino component. 
Assuming this state to be purely active one can estimate both the ^B neutrino flux and 
the Ue survival probability to a much higher precision than currently possible. Finally 
the measured value of the NC rate will have profound implications for the mass and 
mixing parameters of the solar neutrino oscillation solution. 
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A large number of experiments have observed anomalously low solar neutrino 
flux I'^.s.^ compared to the standard solar model (SSM) prediction They are the 
radiochemical experiments on Ga ^ and CI ^ targets as well as the water Cerenkov 
experiments from Super Kamiokande (SK) ^ and Sudbury Neutrino Observatory 
{SNO) SK observes the emitted electron from elastic scattering 

ly + e ^ ly + e, (1) 

while SNO observes it from the charged current process 

Ve + d^p + p + e (2) 

using a heavy water target. Both the experiments probe the high energy tail of 
the solar neutrino spectrum, which is dominated by the neutrino flux. The SK 
elastic scattering process (1) is sensitive to Ve via CC and NC interactions; but it 
also has a limited sensitivity to iy^,r via the NC interaction. On the other hand the 



SNO experiment is looking separately at the NC rate, which has equal sensitivity 
to all the active neutrino flavors, via 

u + d^u + p + n, (3) 

followed by the neutron capture by NaCl. The resulting excited state decays via 
7 emission, which constitutes the NC signal. The SNO experiment is expected to 
report its first NC rate shortly, corresponding to a data sample of ~ 10^ events - 
i.e. similar statistical accuracy as their 1st CC data The purpose of our work is 
to critically analyse how far this data will enhance our understanding of the solar 
neutrino anomaly and its solution. 

By far the most plausible solution of the solar neutrino anomaly is in terms of 
neutrino oscillation, and in particular the oscillation of fg into another active flavor 
Va, which can be any combination of and fr- Indeed we have made considerable 
progress in understanding the solar neutrino anomaly and its oscillation solution 
over the past few months by combining the CC rate from SNO with the SK elastic 
scattering data 6,7,8,9,io pij-stly they disfavour transition into a sterile neutrino 
j^s at the 3(7 level in a model independent way.^ Secondly, assuming only transition 
between Ue and i^a, one can estimate both the flux and the survival probability, 
although with fairly large uncertainties.^ Thirdly one can do a combined analysis of 
the SNO, SK and the earlier radiochemical data assuming SSM and the i^e — * z^o 
oscillation model One sees that the result has narrowed down to only two 

possible solutions, both having large mixing angles - i.e. the so called LMA and 
LOW solutions. We shall see below that each of the above results can be significantly 
sharpened further by including the NC rate from SNO. 

Table 1. The observed solar neutrino rates relative to the SSM predictions are shown along with 
their compositions and threshold energies for different experiments. For the SK experiment the 
Ue contribution to the rate R is shown in parentheses assuming —* Va transition. 



experiment 


R 


composition 


(Mev) 


Ga 


0.584 ± 0.039 


pp(55%),Se(25%),B(10%) 


0.2 


CI 


0.335 ± 0.029 


B(75%),Be(15%) 


0.8 


SK 


0.459 ± 0.017 (0.351 ± 0.017) 


5(100%) 


5.0 


SNO{CC) 


0.347 ± 0.027 


B(100%) 


7.0 



Table 1 lists the observed solar neutrino rates of the Ga, Gl, SK and SNO (CC) 
experiments relative to the corresponding SSM predictions. The compositions of 
the respective solar neutrino fluxes are also indicated along with the threshold 



energies. Assuming the SSM neutrino fluxes and the transition of into an active 
flavor Va one can write the SK elastic scattering rate relative to the SSM prediction 
in terms the survival probability, i.e. 

R'§K=Pee + rP,a=Pee+ril-Pee), (4) 

where r = (7^^ /(^^^^^'^ — ^■^'^ the ratio of v^j^^r to Ue elastic scattering cross- 
sections (|^).^ The resulting value of the Ve survival probability Pee is shown par- 
enthetically in Table 1. The other rates shown in this table are identical to the 
corresponding Pee, since they are sensitive to Ve only. 

As we see from table 1, the SK and SNO experiments are sensitive only to the 
neutrino spectrum. While the shape of this spectrum is predicted with good 
precision by the SSM, there is a large uncertainty in the predicted normalisation, ^ 

= 5.05 X 10*^ (ll;?^) cm-2s-\ (5) 

arising from the uncertainty in the ^Pe -f p — > ®P -I- 7 cross-section. Therefore 
we shall introduce a constant parameter fs to denote the normalisation of the *P 
neutrino flux relative to the SSM prediction. Then the SK elastic scattering and the 
SNO CC and NC scattering rates relative to the corresponding SSM predictions 
are 

RSK = fsPee+fBrPea, (6) 
R-SNO = fsPee, (7) 
PsNO = fB{Pee+Pea), (8) 

which hold for the general case of transition into any combination of Va and 
i^s- It should be noted here that these three measurements do not have identical 
energy ranges. The SK and SNO CC data start from neutrino energies of 5 and 
7 MeV respectively, while the response function of the SNO NC measurement 
extends marginally below 5 MeV. However the SK rate and the resulting survival 
probability show energy independence down to 5 MeV to a very good precision. 
The SNO CC rate shows energy independence as well, although to lesser precision. 
Therefore it is reasonable to assume a common survival probability for all the three 
measurements. 

As recently discussed in, ^ the SK elastic and the SNO CC rates can be com- 
bined to access information on NC scattering, so that in principle there is no new 
information contained in the SNO NC rate. In fact eqs. can be seen to give 

the sum rule 

P-SNO — PsNO + (PsK ^ P^No)/'^^ (9) 

which predicts Rsno ~ l-0±0.24 However the large uncertainty in this prediction 
reflects the low sensitivity of Rg^ to NC scattering, which is weighted by a small 
co-efficient r. On the other hand the expected NC scattering rate from SNO at the 



level of ~ 10^ events should have a similar precision as their present CC rate, i.e. 
±8%, which will ultimately go up to ±5%. Keeping in mind the predicted range of 
^SNO shall assume 

RsNO = 1-0 ± -08, 1.2 ± .08, 0.8 ± .08 (10) 

as three representative values of the NC scattering rate expected from SNO and 
analyse the implications for the solar neutrino anomaly. We shall see below that 
the Rg^o iiiput will lead to qualitative improvement in the results of all the three 
types of analyses mentioned earlier. 
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Fig. 1. The SNO CC and NC rates shown relative to their SSM predictions for the three 
representative values of the latter. The dashed line is the prediction of the pure to Us transition 

General Analysis of i/e Va,s Transitions: Let us start with a model indepen- 
dent analysis of transition into an active neutrino flavour v^- One sees from eqs. 



that the observed excess of R'§x '^^^^ t^"^ ^sno Table 1 constitutes a 3cr 
signal for i^e ~> i^a transition or equivalently a 3cr signal against a pure I'e i^s 
transition Similarly an observed excess of Rg^p ® '^^^^ ^sno (0) ^^^^ constitute 
a model independent signal for the ~> i^a transition P^a- Fig. 1 compares the 
^SNO values of eq. ( p^ ) with the current value of Rs^q- shows that an observed 
value of Rg^o range of 0.8 — 1.2(±.08) will constitute of 5 — 10 ct signal for 

~^ transition i.e. transition of Ve into a state of atleast partly active neutrino. 




Fig. 2. Best fit value of tfie °B neutrino flux shown along with the Icr and 2a limits against the 
model parameter sin^ a, representing Ue transition into a mixed state {uasina + i^acosa). The 
dashed line denote the ±2(r limits of the SSM. 



Unfortunately one cannot get a model independent estimate of the Ve — > i^a 
transition probability Pea, since the above excess corresponds to the product fsPea- 
In fact it is evident from eqs. (^|-^) that these two quantities can not be separated 
with or without the Pg^Q input. In other words for the general case of v,, transition 
into a mixture of active and sterile neutrinos, 



Va sin a + cos a, 



(11) 




Fig. 3. Same as Fig. 2 but including the NC rate to the fit along with the SK elastic and SNO 
CC rates. The dashed lines denote the ±2(7 limits of the SSM . 

one can rewrite eqs. (||) and ^ as 

RSK = fBPee+ f BT sin'' a{l-P,e), (12) 
RSNO = fsPee+fBSm^ail-Pee). (13) 

Then it will not be possible to separately estimate the parameters Jb and sin^ a. 
One can only determine the following combination in two different ways, i.e. 

sin^ a{fB - R'^^o) = iRlK-RsNc)^ (14) 
sin' aijB - R'^So) = RsSo-RsSo- (15) 

A two parameter fit to (|lj) was done in ref.® to determine the allowed contour in 
the fs — sin^ a plane. We shall instead treat sin^ a as a model parameter. And 
for each input value of sin^ a we shall determine fg first from eq. ( |l^ and then 
from a weighted average of eqs. (|lj) and (p^). The results are shown as lines 
corresponding to the central value of fs along with the la and 2a boundaries in 



the Jb — sin^ a plane in Figs. 2 and 3. The advantage of this procedure is that 
for sin^ a = 1 the la and 2(t ranges of /s shown correspond to those of the pure 
i^e i^a transition model without any sterile neutrino, which would not be the case 
for the two parameter fit 

A comparison of Figs. 2 and 3 shows the enormous improvement in precision 
one expects by including the NC rate from SNO. The vertical lines indicate the 
2(7 limits on Jb from the SSM. Combining the two limits one could rule out models 
having sin^ a < 0.3 (i.e. singlet neutrino component cos^ a > 0.7) at 2a level for 
R^c^ = 0.8 ± .08, while R^^o = 1-2 ± -08 would rule out all those singlet neutrino 
models with sin^ a < 0.6 (cos^ a > 0.4). However within the present uncertainty 
of the SSM value of the neutrino flux it is unlikely to rule out models with a 
i/fi component cos^ a < 0.3. Nor is it likely to set any upper limit on sin^ a, which 
would show incompatibility with pure Ve — > i^a transition. 

Analysis of pure i/^ — > i^a Transition: Assuming no sterile neutrino the eqs. (^) 
and (^) reduce to 

RSK = fBPee + fBr{l-Pee), (16) 

RsNO = Ib- (17) 

In this case it is possible to estimate both /b and Pee from the SK elastic and 
the SNO CC rates even without the NC rate from SNO J Fig. 4 shows the 1 
and 2 sigma contours of this solution along with the corresponding ones obtained 
by including the Rg^o input. Our result in the former case is in good agreement 
with that of ref. ^, although the two methods are not identical. Thanks to the 
effective energy independence of the survival probability the result is insensitive 
to the difference between the energy thresholds of the SK elastic and the SNO 
CC rates. Unfortunately the resulting /b has large error, which reflects the low 
sensitivity of i?^^ to /b (eq. ^6|). Interestingly both the central value and the error 
bar of Jb are practically the same as those of the SSM (shown for comparison near 
the right scale) . The error in /b propagates into Pee due to a strong anticorrelation 
between the two quantities as their product is well constrained by Rg^Q Fig. 4 
shows that the inclusion of the NC rate from SNO will result in a reduction in the 
uncertainty of the neutrino flux /b and the corresponding survival probability 
Pee by about a factor of 2.5. 

The i^e Va oscillation solutions to the Solar Neutrino Anomaly: Finally 
we shall fit the global solar neutrino data with a standard two-family (i^e Va) 
oscillation model assuming the SSM fluxes, but with one difference. Instead of 
Rgj^ and Rsno shall fit the ratios Rsk/P-sno ^^'^ P-sno/^SNo- ^^e 
from eqs. (|7l |l^ , [l7| ) the ^B neutrino flux Jb factors out from these two ratios. Thus 
the result becomes immune to the large uncertainty in the SSM value of the ^B 
neutrino flux (||). We shall include the 19 + 19 day-night spectral points from SK, 
but with free normalisation to avoid double counting J'^'^^ We shall also include 




Fig. 4. The 1 and 2(t contours of solutions to the neutrino flux fs and the Ue survival 
probability Pee assuming i/e to i/a transition. The size of the SSM error bar for fg is indicated 
on the right. 



the combined Ga rate of Tabic 1 in the fit. However we shall exclude the CI rate, 
since the experiment has not been independently calibrated. Besides the apparent 
rise of the CC rate between the CI and SK/SNO energies is in conflict with the 
LMA and LOW solutions, which are strongly favoured by the global fit.^*' Fig. 5 
shows the 90%, 95%, 99% and 99.73% (3cr) CL contours of the fits for the three 
representative values of R^no (0)' which is the same as /b (eq. We have 

checked that all the three fits have equally good x^in (~ 30) and goodness of fit 
(~ 80%). 

The left panel of Fig. 5 represents Rf/cF = Ib = 0.8 ± .08. It corresponds to 
enhancing the survival probability Pee at SK/SNO from 0.35 to nearly 0.45, while 
Pee at Ca energy goes up by only 2%. Thus the energy dependence of Pee between 
the Ga and SK/ SNO energies become very mild. Consequently the solution favours 
large mixing angle, going upto maximal mixing, where it remains valid over two large 
mass ranges around Sm'^ = 10"'' eV^ and lO"'^ eV^ These are the so called LM A 
and LOW solutions. In the present context however it will be more appropriate 
to call them HIGH and LOW solutions, since both of them correspond to large 
mixing angles. The middle panel represents Rg^o = /b = 1-0 ^ -08, which means 
that the survival probability Pee at Ga and SK/SNO energies corresponds to those 
shown in Table 1. Since such a large energy dependence cannot be explained by the 
earth matter effect, the LOW solution is only allowed at the 3cr level. On the other 
hand the HIGH solution occurs at the upper edge of the A4SW region in SnP' /E. 
Thus the survival probability goes down from Pee — 1 ~ 5 sin^ 29 > 0.5 above 
the MSW region to Pee — sin^ 9 inside it as one moves up from Ga to SK/SNO 
energies .^^ The right panel corresponds to Rg^o = Jb = 1-2 ± .08, which means 
that the survival probability at SK/SNO energies goes down further by a factor 
of 1.2. Consequently the LOW solution disappears completely while the HIGH 
solution moves to a lower mixing angle. Thus the measured value of the NC rate at 
SNO can have a profound effect on the mass and mixing angle of the solar neutrino 
oscillation solution. 

Summary: 

In anticipation of the first NC rate from SNO we have analysed how this data 
will sharpen our understanding of the solar neutrino anomaly. For a quantitative 
analysis we have chosen three representative value of this rate, Rs^o ~ 0.8, 1.0 and 
1.2(±.08). They span the ±lcr range of this quantity as estimated from the SK 
elastic and SNO CC rates. The main results are listed below. 

(i) It will provide a 5 — lOcr signal for Ve transition into an active flavor (or 
against a pure Ve transition into a sterile neutrino i>s)- 

(ii) However for transition into a mixture of Va and Vg, we need to know the 
neutrino flux to constrain the size of the sterile component. If we assume the 
SSM prediction for this flux then R^^q = 1.2 ± .08 (0.8 ± .08) would imply 
the sterile component to be < 40% (70%) at the 2a level. 
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Fig. 5. The Ve —* Va oscillation solutions to the Ga rate, SK day-night energy spectra along with 
the SK and SNO (CC) rates, both normalised by the SNO (NC) rate. 



(iii) Assuming a pure Ve ^ transition one can combine Rg^p "^i^h Rsno ^^'^ 
^'sK to determine both the neutrino flux and the v,. survival probability 
to much higher precision than is possible now from the latter two data. 

(iv) Finally one can do a i^e ^ fa oscillation model fit to the global solar neutrino 
data assuming the SSM fluxes, but replacing the i?^^ ^sno ^^'^ ratios 
R'sk/ ^SNO ^'^^ ^SNo/-^SNO^ which are independent of the neutrino flux. 
For R^gfj = 0.8 ± .08 we get both the LMA (HIGH) and LOW solutions 
covering large ranges of mass and mixing angle, including maximal mixing. 
On the other hand Rg^Q > 1 strongly disfavours the LOW solution, while 
the HIGH solution is restricted to a small patch in mass and mixing angle 
excluding maximal mixing. Thus the measured value of Rg^o '^^^ have a 
profound effect on the mass and mixing parameters of the solar neutrino 
oscillation. 
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